Observation of Multi-Gap Superconductivity in GdO(F)FeAs by Andreev 

Spectroscopy 
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We have studied current-voltage characteristics of Andreev contacts in polycrystalline 
GdOo.88Fo.i2FeAs samples with bulk critical temperature T c = (52.5 ± 1) K using break-junctions 
technique. The data obtained cannot be described within the single-gap approach and suggests the 
existence of a multi-gap superconductivity in this compound. The large and small superconducting 
gap values estimated at T = 4.2K are Al — 10.5 ± 2 meV and As = 2.3 ± 0.4 meV, respectively. 

PACS numbers: 74.70.Xa, 74.20.Fg, 74.25.-q, 74.45. +c, 74.62.Dh 



Novel superconducting compounds of 1111 family 
based on rare-earth oxypnictides REOFcAs (RE = La, 
Sm, Gd etc.) [II, [3 are currently in the focus of research 
interest. Some of their features such as layered structure 
and spatial separation of the carrier reservoir layers and 
the superconducting pairing layers are similar to those of 
cupratcs. However, many other properties differ substan- 
tially and promise new interesting physics At present, 
the key issues under investigations are the effect of vari- 
ous types of doping, pairing mechanism, symmetry of the 
order parameter, quasiparticle energy spectrum, and the 
superconducting energy gap(s). 

The stoichiometric compounds of the 1111-family are 
antiferromagnetic metals with spin density wave ground 
state [i| . Partial deficiency of oxygen or fluorine substi- 
tution for oxygen induces superconductivity in the FeAs- 
layers. Replacement of rare-earth elements also affects 
the superconducting critical temperature, T c . In par- 
ticular, T c of Gd-bascd oxypnictide may be lowered by 
partial replacement of Y for Gd [f| or gained up by 
introducing Th instead of Gd @. T c =56 K found in 
Gdo.8Tho.20FcAs compound is today the highest one for 
iron-based superconductors. 

According to band structure calculations, the total 
density of states at the Fermi level N(0) is formed mainly 
by Fe 3d-states 0-01 ■ As shown in Ref. 0, the T c values 
for different iron-based superconductors correlate with 
N(0), thus giving support to the BCS-like coupling in 
these compounds. 

The theoretically calculated Fermi surface for 1111- 
system (ll| - [l3| consists of quasi-two-dimensional (2D) 
hole sheets centered at the T point and two electron 
sheets at the M points of the first Brillouin zone. Within 
the so called minimal two-band model, these four bands 
may be considered as two effective 2D bands [14, 



Correspondingly, many of the available theoretical and 
experimental data indicate that iron-based layered ma- 
terials arc multiband superconductors with s-type sym- 
metry of the order parameter 0). Knight shift measure- 
ments in 1111-class compounds jiol ] have proven unam- 
biguously the spin-singlet type pairing in these materials. 
Several data were reported in favor of s ± [13, E3 or s ++ 
order parameter symmetry (l9l | , making the experimental 
situation regarding 1111-compounds uncertain. 

The magnitude and structure of the superconducting 
gap A is intimately related to the pairing mechanism. 
ARPES measurements are not sensitive enough to resolve 
unambiguously such fine details as A, on the scale of a 
few meV, that makes this parameter accessible nearly 
exclusively from point contact spectroscopy, such as 
scanning tunneling spectroscopy (STS), tunneling- and 
point-contact Andreev reflection (PCAR) spectroscopy 
(the latter in the regime of SN-, or symmetrical SNS- 
junctions). The available experimental reports are how- 
ever rather inconsistent for 1111 class compounds [20I — 
HH, even for the most intensively studied SmO(F)FeAs. 
Various types of conclusions have been reported includ- 
ing d-wavc like, single gap-like, and multi-gap behavior. 

The ambiguity of the experimental information is 
partly due to an inevitable inhomogeneity of the 1111- 
type polycrystalline samples and lack of large size 1111- 
type single-crystals suitable for these measurements. An- 
other cause for the divergency of the point contact spec- 
troscopy data is inherent in those experimental tech- 
niques, where the sample surface is not cleaved in high 
vacuum or cryogenic environment. In order to resolve 
the experimental ambiguity, evidently, novel sets of com- 
prehensive experimental data are needed, which would 
comprise self-consistency check, substantial statistics and 
provide local probing at various points of the in-situ 
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cleaved surfaces. 

Here we report the superconducting gap measurements 
in nearly optimally doped GdOo.8sFo.i2FeAs samples by 
SNS Andreev spectroscopy using the break-junction tech- 
nique ■ Until now these measurements have not been 
done for Gd-1111, an analogue to Sm-1111 with approxi- 
mately the same T c ~ 53K. The break junction tech- 
nique opens a nice opportunity to prepare in helium 
atmosphere, at liquid 4 He temperatures, clean surfaces 
forming Andreev contact. Another advantage is a possi- 
bility of fine mechanical readjusting the contact during 
experiment, that enables to collect multiple data from 
different local areas of the same sample. Using this 
method we have unambiguously detected the presence 
of two superconducting gaps, whose best fit values aver- 
aged over about 30 spectra are A^ = 10.5 ± 2meV and 
A s = 2.3 ± 0.4 meV at T = 4.2K. 

Polycrystalline samples GdO(F)FeAs were prepared by 
high pressure synthesis (2^|. The chips of high purity 
Fe, and powders of single-phase FeF 3 , Fe203, and GdAs 
were mixed together in the nominal ratio and pressed 
into pellets of 3mm diameter and 3mm height. The pel- 
lets were placed in boron nitride crucible and synthe- 
sized at pressure of 50 kb and temperature 1350°C during 
60 min. The X-ray diffraction pattern averaged over the 
sample area showed a polycrystalline compound with a 
dominating desired 1111-phase (with lattice parameters 
a = 3.902(2)A, c = 8.414(5)A ) and an admixture of in- 
cidental FeAs and Gd2C>3 phases. The subsequent local 
EDS analysis (JSM-7001FA) has revealed that the inci- 
dental phases are concentrated in grains of about lmkm 
size which are scattered in the bulk majority phase. This 
fact opens a possibility to probe properties of the true 
majority phase using local techniques, such as PCAR. 
Superconducting properties of the samples were tested by 
measurements of temperature dependence of ac-magnetic 
susceptibility and resistivity R(T). Both showed a sharp 
superconducting transition in our polycrystalline samples 
with T c rj 52. 5K (the latter value was defined at a max- 
imum of dR(T) /dT-cmve). Figure 1 shows typical tem- 
perature dependence of resistance and its derivative. 

For point-contact spectroscopy we used two meth- 
ods: (i) multiple Andreev reflections spectroscopy of 
individual superconductor-constriction-superconductor 
Sharvin-type contacts (26l - [28| and (ii) intrinsic Andreev 
specrtoscopy of stack contacts that usually exist due to 
the presence of steps and terraces on clean cryogenic 
cleaves in layered crystals. 

Thin plates of about 2 x 1 x 0.12 mm 3 in size were 
cut from the synthesized pellets. At room temperature, 
the plate-like sample was mounted onto an elastic bronze 
holder and the two current and two potential leads were 
attached to the sample by liquid In-Ga alloy. The holder 
with the sample was placed in the measuring cell and 
cooled down to 4.2 K. A microcrack in the sample was 
generated by precise bending the sample holder at T = 




FIG. 1: Fig. 1. Superconducting transition for polycrystalline 
GdOo.88Fo.i2FeAs sample measured prior a microcrack for- 
mation (dots). The bulk T c = (52.5 ±1) K was determined at 
a maximum in dR(T) / dT-curve (solid line). 



4.2K using a micrometric screw. 

Current- voltage dependence, I{V), and its derivative, 
dI(V)/dV, were measured automatically using the 16-bit 
AT-MIO-16X (National Instruments) digital board. The 
amplitude of a low-level 820 Hz modulation voltage at 
potential leads of a sample was maintained stable using 
a lock-in nanovoltmctcr (operated as null- detector) and 
a computer controlled digital bridge with a proportional- 
integral-derivativc feedback signal. As a result, the dif- 
ferential conductance of a contact was proportional to 
the amplitude of the ac feedback current through the 
contact. 

Figure 2 represents I(V), dI{V)/dV and d 2 I{V)/dV 2 
characteristics for individual Andreev (SNS) break- 
junction in polycrystalline GdOo.8sFo.i2FeAs sample 
measured at T = 4.2K. The observed experimental IV- 
curves are typical for the clean classical SNS-contacts 
with excess-current 



characteristics [27|, 12911, therefore, 
the theoretical model of Kummcl et al. [27| is sup- 
posed to be applicable to our break-junctions. Accord- 
ing to the Kummel model, the TV-characteristics at low 
bias voltages should show a subharmonic gap structure 
(SGS) with a series of dips in the dynamic conductance 
dI(V)/dV at bias voltages 



V n = 2A/en, 



(1) 



with an integer n = 1,2..., due to multiple AR effect. 
For a two-gap superconductor, two independent SGSs 
corresponding to the large A^ and small As gaps are 
anticipated. 
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They may be caused by excitation of collective modes 
and require additional studies. 
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FIG. 2: Fig. 2. /(V), dI{V)/dV and d 2 7(y)/dV 2 -curves for 
a single SNS-contact 1D06 at T = 4.2K. Background (a poly- 
nomial function) is subtracted. The set of dips in the dif- 
ferential conductance at bias voltages V„l = 2Ai/en (verti- 
cal solid lines) determines the energy of the large supercon- 
ducting gap, Al ~ llmeV. Peculiarities on the dI(V)/dV 
and d 2 I{V)/dV 2 -curves, marked by dashed lines, indicate the 
presence of small superconducting gap, As « 2.5 meV. 
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FIG. 3: Fig. 3. Differential conductance of the SNS contact 
3D0I in GdOo.88F .i 2 FeAs sample at T = 4.2K. Background 
is subtracted. Two differential conductance dips define the 
energy of the small superconducting gap As ~ 2.15 meV. 
The anticipated bias voltages V n s = 2 As /en are depicted by 
vertical dashed lines. 



The dips labeled on Figure 2 as Ul = 1 and 2 on SGS 
reflect the large gap; they are marked with vertical solid 
lines. The singularities shown by vertical dashed lines 
cannot be attributed to the large gap and, therefore, may 
reflect the existence of a small gap As « 2.5 meV. Com- 
paring the result for the large gap (Fig. 2) with Eq. (1) 
one can easily obtain Ax, ~ llmeV. 

By readjusting the contact, we could observe clear sets 
of dips on dI(V)/dV curve due to either large- or small- 
gaps, or even to both (as shown in Fig. 2). Figure 3 
expands the small bias range, where the subharmonics of 
the small gap are clearly seen in dl /dV curves for a single 
SNS-contact EL3D01. For clarity, the smoothly varying 
background is subtracted. Using expression (1), from 
the set of differential conductance dips on the dI(V)/dV 
curve we determined the energy of the small supercon- 
ducting gap As = 2.15 meV. In Fig. 3, one can also see 
two extra features at w 3mV and 3.7 mV, which can not 
be attributed to either large or small gap subharmonics. 



The peculiarities in the dI(V)/dV and d 2 I(V)/dV 2 - 
characteristics shown in Figs. 2 and 3 clearly mani- 
fest the existence of two gaps in our GdO(F)FcAs sam- 
ples. This behavior is similar to that observed earlier 
in the multi-band superconductor Mgi_ :c AL J ;B2 and 
in LaOo.gFo.iFeAs [31| (an analog to our Gd-1111 sam- 
ple with somewhat lower T c w 28K). The sharpest SGS 
(like those in Figure 3) may be usually observed only on 
dI(V) I dV-characteristics of Andreev contacts of the high 
quality and of small size, comparable to the quasiparti- 
cles mean free path (ballistic limit) {28|. For such a case, 
a number of the observable gap peculiarities (up to 4 in 
some samples) facilitates interpretation of the multigap 
subharmonic structure. 

Our experimental data are summarized in Figure 4 
where the normalized to a single junction bias voltages 
VnL,S for four microcontacts are plotted versus l/n^s. 
According to expression (1), such dependences have to 
fall onto straight lines passing through zero. This is in- 
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FIG. 4: Fig. 4. Normalized bias voltages = 2AL,s/en ver- 
sus l/riL,s for the studied SNS-arrays. The averaged values 
of the superconducting gaps are Al — (10.5 ± 2)meV and 
As = (2.3 ± 0.4)meV. Solid lines are guides to the eye. 



deed fulfilled for in all samples; for As, this however, 
could be verified only in sample EL3D01, because a rich 
picture of reproducible features detected at low bias volt- 
ages impeded the analysis for other samples. 

Based on the data obtained we conclude on the ex- 
istence of two distinct superconducting gaps with ener- 
gies Al = (10.5 ± 2) meV and A s = (2.3 ± 0.4)meV at 
T = 4.2K in GdOi-^F^FeAs sample. The reproducibil- 
ity of two SGSs detected at d/(l^)/dl^-characteristics of 
various Andrccv arrays, formed by the break-junction 
technique support this conclusion. In some cases we ob- 
served extra features in dl/dV- curves signalling the ex- 
istence of the 3rd, smaller gap, Ass ~ 1 meV. 

Using the determined gap energies and bulk T c = 
(52.5 ± 1)K, one can estimate 2A/fcsT c ratio. For the 
large gap, our experimental data lead to 2Al//sbT c = 
(4.8 ± 1.0) that exceeds the standard BCS value, 3.52, for 
single-gap superconductors in the weak coupling limit. 
This fact together with rather conventional exponent 
value for Fe isotope effect [32| resembles the BCS - 
model behavior with strong electron-phonon coupling. 
At the same time, the 2A/ksT c ratio for the small gap 



2As/ksT c ss 1.1 <C 3.52 suggests that the "weak" super- 
conductivity may be induced by intcrband coupling, due 
to fc-space internal proximity effect between two conden- 
sates, where the large gap condensate plays the "driving" 
role. In particular, similar situation is believed to be re- 
alized in MgB 2 Hi and LaO .gF .iFeAs [H|. 



TABLE I: Table 1. Summary of the A values (in meV) 
measured for 1111-family REO(F)FeAs compounds by point 
contact Andreev reflection (PCAR), break-junction Andreev 
reflection (BJ), scanning tunneling spectroscopy technique 
(STS), and ARPES 



RE 


Tc(K) 


method 


2A L 


2A S 


ref. 


Gd 


53 


BJ 


10.5±2 


2.3±0.4 


this work 


Sm 


53 


ARPES 


15±1.5 


no 


[34] 


Sm 


52 


PCAR 


18±3 


6.15±0.45 


[35] 


Sm 


52 


PCAR 


19 


5.7 


[36] 


Sm 


52 


STS 


8-8.5 


no 


[21] 


Sm 


51.5 


PCAR 


20 


6.6 


[36] 


Sm 


51 


PCAR 


10 


4 


[22] 


Nd 


51 


PCAR 


14±1 


6 ±1 


[37] 


Nd 


51 


PCAR 


12.5±0.5 


6.3±0.3 


[38] 


Nd 


48 


BJ,STS 


7-10 


no 


[33] 


Tb 


45 


PCAR 


8.8 


5 


[39] 


Nd 


45 


PCAR 


11±2 


5 ±1 


[40] 


Sm 


42 


PCAR 


15±1 


4.9±0.5 


[35] 


Sm 


42 


PCAR 


6.7±0.1 


no 


[18] 



The presence of the large superconducting gap char- 
acterized by lAj^jkBTc > 3.52 in the 1111-family com- 
pounds REOFcAs (RE = La, Sm, Nd) was confirmed 
by tunneling spectroscopy using break-junction tech- 
ni que 1331 . point-contact Andreev reflection spectros copy 
pel l22l I35I4I0I ] . scanning tunneling spectroscopy pH l33f. 
and angle-resolved photocmission spectroscopy (ARPES) 
34j (see Table 1). To the best of our knowledge, there 
is no other available data for Gd-1111. Therefore, we 
compare in Table 1 our data for Gd-1111 with other data 
available for Sm-, Nd- and Tb-1111 superconductors with 
similar T c . We emphasize rather good agreement between 
2Ai/fcsT c values determined from our study and those 
from STS 21], break-junction measurements 
some PCAR- measurements fH&Lil. 



and 



As to the small gap, there is evidently a sizeable spread 
in its value (2 — 5)mcV observed in different experiments. 
We can not also exclude that the spread may be caused 
by the existence of three- gap superconductivity in 1111- 
system with multiple-sheet Fermi surface [l3 , 41 1 ; for the 
smallest gap, we estimate 2A$ /ksTc m 0.5. 

In conclusion, we have studied the I(V)- and 
d/(F)/(il / -characteristics at T = 4.2K for vari- 
ous SNS Andreev break-junctions in polycrystalline 
GdOo.88Fo.i2FeAs samples with bulk critical tempera- 
tures T c m 52. 5K. The obtained characteristics do not fol- 
low the standard single-gap model behavior. Two clearly 
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observed independent subharmonic gap structures point 
at the existence of two distinct superconducting gaps, 
A L = (10.5 ± 2)meV and A s = (2.3 ± 0.4)meV deter- 
mined at T = 4.2K. The estimated 2Ai J /ksT c ratio ex- 
ceeds the standard BCS value, 3.52, for single-gap super- 
conductors and weak-coupling limit while for the small 
gap the 2A/T c -ratio 2A s /k B T c <C 3.52. 
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by the Russian Foundation for Basic Research and by the 
Russian Ministry for Education and Science. 
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